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Through the introduction of a class of appropriate translationally invariant trial wave functions, 
we show that the strong correlations in the lowest Landau level (LLL) reflect in finite systems the 
emergence of intrinsic point-group symmetries associated with rotations and vibrations of molecules 
formed through particle localization. This molecular description is universal, being valid for both 
bosons and fermions, for both the yrast and excited states of the LLL spectra, and for both low and 
high angular momenta. This physical picture is fundamentally different from the "quantum-fluid" 
one associated with Jastrow-type trial functions. 
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Motivation. — Following the discovery [l[ of the frac- 
tional quantum Hall effect (FQHE) in two-dimensional 
(2D) semiconductor heterostructures under high mag- 
netic fields (B), the description of strongly correlated 
electrons in the lowest Landau level (LLL) developed into 
a major branch of theoretical condensed matter physics 
@, ! i, B i 0, i i M EH, H M H 0- Early 
on, it was realized that the essential many-body physics 
in the LLL could be captured through trial wave func- 
tions. Prominent examples are the Jastrow-type Laugh- 
lin 0], composite fermion [f|, and Moore and Read's 
7] Pfaffian functions, representing quantum-liquid states 
2j. More recently, the field of semiconductor quantum 
dots [HI helped to focus attention on finite systems with 
a small number (N) of electrons. Theoretical investi- 
gations of such finite systems led to the introduction of 
"crystalline" -type LLL trial functions referred to as ro- 
tating electron molecules (REMs) (l2. 14\. 



Most significantly, the burgeoning field of trapped ul- 
tracold neutral gases has generated recently an unparal- 
leled interest regarding the fundamental aspects (includ- 
ing the appropriateness of trial functions) of stron gly cor- 
related states in the lowest Landau level [l7], HRIla. I20L 
HI, HI 0, 0, HI, [11; in this case, the LLL manifold 



of degenerate single-particle orbitals arises as a result of 
the rapid rotation (with rotational frequency Q) of the 
trap. Furthermore, it is anticipated that small assemblies 
of ultracold bosonic atoms will become technically avail- 
able in the near future 12 711 and that they will provide 
an excellent vehicle [22I l24l Eol [27} for experimentally 
reaching exotic physical behavior beyond the mean-field 
Gross-Pitaevskii regime and for testing the rich variety 
of proposed LLL trial wave functions. 

A universal description of the full LLL spectra (in- 
cluding both yrast [28| and all excited states), how- 
ever, is still missing. To remedy this, a unified the- 
ory for the LLL spectra of a small number of parti- 
cles valid for both statistics (i.e., for both bosons and 
fermions) is introduced in this paper. The LLL spec- 
tra are shown to be associated with fully quantal ro- 



vibrational molecular (RVM) states, i.e., states described 
by trial wave functions akin to the REM functions of Ref. 
[12j . It is remarkable that the numerical results of the 
present theory agree within machine precision with exact- 
diagonalization (EXD) results, including energies, wave 
functions, and overlaps. This behavior outperforms the 
behavior of all other trial functionSjincluding those of 



the composite- fermion view ]1| 1|| 2t| . The RVM func- 



tions and corresponding point-group symmetries express 
the emergent many-body intrinsic structure of the highly 
correlated LLL states. 

Theory. — The RVM functions have the general form 
(within a normalization constant): 



RXM 



(n l ,n 2 )Q%\0>, 



(1) 



where (714,712) indicates the molecular configuration 
(here we consider two concentric rings) of point-like 
particles with n\ (77,2) particles in the first (second) 
ring. The particles on each ring form regular polygons. 
The index RXM stands for either REM, i.e., a rotat- 
ing electron molecule, or RBM, i.e., a rotating boson 
molecule. $£ XM (77,1, 77,2) alone describes pure molecu- 
lar rotations associated with magic angular momenta 
C = Cq + riiki + 7i2fc2, with k%, k 2 being nonnegative 
integers; Co = N(N — l)/2 for electrons and Co = 
for bosons. The product in Eq. (fT]) combines rotations 
with vibrational excitations, the latter being denoted by 
Q™, with A being an angular momentum; the superscript 
denotes raising to a power m. Both $2 XM and are 
homogeneous polynomials of the complex particle coordi- 
nates z\, Z2, ■ ■ ■ , zn, of order C and Am, respectively. The 
total angular momentum L = C+Xm. Q™ is always sym- 
metric in these variables; <J)^ XM is antisymmetric (sym- 
metric) for fermions (bosons). |0 >= YiiLi ex P[ — z i z i*/2]; 
this product of Gaussians will be omitted henceforth. 
The analytic expressions for the $^ EM (for fully po- 



larized electrons) were derived in Ref. [12j employing a 
two-step method: (i) First a single Slater determinant 
[that breaks the rotational (circular) symmetry] was con- 
structed using displaced Gaussians as electronic orbitals, 



2 



symmetric polynomials: 



u(z, Zj) = —= cxp[-|z - Zj\ 2 /2] exp[-i(xYj + yXj)].(2) 
V 71 " 

The phase factor is due to the gauge invariance. z = 
x — iy, and all lengths are in dimensionless units of 
IbV%, w ith the magnetic length Ib = yjh/ (m e oj c ); 
uj c = eB/(m e c) is the cyclotron frequency. The cen- 
ters Zj = Xj + iYj, j = 1, 2, . . . , N of the Gaussians are 
the vertices of the regular polygons in the (711,71.2) geo- 
metric arrangement, (ii) A subsequent step of symmetry 
restoration was performed using the projection operator 
V{C) = ± dje^-V, where L = k is the to- 

tal angular momentum operator; this yielded trial wave 
functions with good total angular momenta C [121 Il4j . 

Analytic expressions for the $^ BM (for spinless bosons) 
can also be derived using the two-step method. Natu- 
rally, in the first step one constructs a permanent out 
of the orbitals of Eq. ((2]); one also uses the equivalence 
lu c — > 2Q between the cyclotron fre quen cy (electrons) and 
the rotational frequency (bosons) [14|]. The expressions 
for any number N of bosons and any molecular config- 
uration (rii, 712, . . . , n q ) will be pressented in Ref. [301 ] . 
Here we present as an illustrative example the simpler 
case of N = 3 bosons having a (0, 3) one-ring molecular 
configuration. One has (within a normalization constant) 



~h J31 

7 2 > *3 J i 



<I>* BM (0,3) = £ C(h,l 2 ,l 3 ) Perm^ 1 
o<;i<i2<i 3 

(3) 

where the symbol "Perm" denotes a permanent with el- 
ements z\ 3 , i,j = 1,2,3; only the diagonal elements are 
shown in Eq. ((3|). The coefficients were found to be: 



,1<1<J'<3 



27T(7 i -Z J ) 



(4) 



where 1 < M < 3 denotes the number of different indices 
in the triad (/1, l 2 , 13) and the p^s are the multiplicities of 
each one of the different indices. For example, for (1,1,4), 
one has M = 2 and p\ — 2, p 2 = 1. 

The $R EM expressions for electrons in a (0, N) or a 
(1, N — 1) configuration are given by Eqs. (2) and (4) of 
Ref. [HI, respectively. For electrons (1) M = N in all 
instances and (2) a product of sine terms replaces the 
sum of cosine terms appearing in Eq. ([?]). 

We note that ^ XM (n 1 ,n 2 ) — for both bosons and 
electrons when C ^ Co + niki +n.2&2- This selection rule 
follows directly from the point group symmetries of the 
(711,712) molecular configurations. 

The vibrational excitations Q\ are given by the same 
expression for both bosons and electrons, namely, by the 



N 

Qx = E( z » 

i=l 



(5) 



where z c — (1/N) J2i=i z * ^ s the coordinate of the center 
of mass and A > 1 is a prime number. Vibrational exci- 



tations of a similar form, i.e., Q\ = 



N 



(and certain 



other variations), have been used earlier to approximate 
part of the LLL spectra. Such earlier endeavors pro- 
vided valuable insights, but overall they remained incon- 
clusive; for electrons over the maximum density droplet 
[with magic C — Co], see Refs. [§] and Q; for electrons 
over the v = 1/3 (y = Co/ X) Jastrow-Laughlin trial func- 
tion [with magic C — 3Cq], see Re f. II Oil : and for bosons 
in the range < L < N, see Refs. [ll2a[21[. 

The advantage of expression (JSJ 31] is that it is 
translationally invariant (TI), a property also shared 
by both the $g BM and $g EM trial functions. In the 
following, we will discuss illustrative cases, which will 
demonstrate that the molecular trial functions of Eq. 
([l} provide a correlated basis that spans the TI sub- 
space 0, [20I . l32| of nonspurious states in the LLL 
spectra. The dimension D T1 (L) of the TI subspace 
is much smaller than the dimension D EXD (L) of the 
exact-diagonalization (EXD) space spanned by the un- 
correlated determinants Det[zJ x , . . . , z^] or permanents 
Perm^ 1 , . . . ,z l £]. The remaining D EXD (L) - D TI (L) 
states are spurious center-of-mass excitations (generated 
by applying <5™) whose energies coincide with those ap- 
pearing at all the other smaller angular momenta [a]. 
Thus D TI (L) = D EXD (L) - D EXD (L - 1); see Tables U 
andHU 

Three spinless bosons — Only the (0, 3) molecular con- 
figuration and the dipolar A = 2 vibrations are at play 
(as checked numerically), i.e., the full TI spectra at any 
L are spanned by the wave functions 



RBM 
3fc 



(0,3)Q 2 "^{fc,m}, 



(6) 



with k : m = 0,1,2,..., and L = 3fc + 2m; these states are 
always orthogonal. This represents a remarkable analogy 
with the case of N = 3 electrons (see below). 

Table Q] provides the systematics of the molecular de- 
scription for the beginning (0 < L < 12) of the LLL 
spectrum. There are several cases when the TI subspace 
has dimension one and the exact solution coincides with 
a single {k,m} state. For L = the exact solution co- 
incides with $^ BM = 1 (QO = 1); this is the only case 
when an LLL state has a Gross-Pitaevskii form, i.e., it 
is a single permanent [see |0) in Eq. {!])]. For L — 2, we 
found $^ act oc Q 2 (for the subscript [i], see caption of 
Table HI, and since [see Eq. ©J Q2 oc {z\ — z c )(z 2 — z c ) + 
(zi — z r ) (Z3 — z r ) + (Z2 — Zr) (zz — z r ) , this result agrees with 
the findings of Refs. [2^, 34 1 concerning ground states of 
bosons in the range < L < N. For L = 3, one finds 
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TABLE I: Spectra of three spinless bosons. Second column: Dimensions of the EXD and the nonspurious TI (in parenthesis) 
spaces. Fourth to sixth columns: Matrix elements [in units of g/(ivA 2 ), A = y/h/(mQ)] of the repulsive contact interaction 
gS(zi-Zj) between the nonspurious states {k, m} [see Eq. (J6j ] . The total angular momentum L = 3k + 2m. Last three columns: 
Total energy eigenvalues from the diagonalization of the associated matrix of dimension D TI (L). There is no nonspurious state 
with L = 1. The full spectrum at a given L is constructed by including, in addition to the listed TI total-energy eigenvalues 
[D TI (L) in number], all the energies associated with angular momenta smaller than L. An integer in square brackets indicates 
the energy ordering in the full spectrum, with [1] denoting an yrast state. Seven decimal digits are displayed, but the total 
energies agree with the EXD ones within machine precision. 



L 


D EXD (D TI ) 


{k, m} 


Matrix elements 






Total energy 


eigenvalues (TI) 





1(1) 


{0,0} 


1.5000000 






1.5000000[1] 




2 


2(1) 


{0,1} 


0.7500000 






0.7500000[1] 




3 


3(1) 


{1,0} 


0.3750000 






0.3750000[1] 




4 


4(1) 


{0,2} 


0.5625000 






0.5625000[2] 




5 


5(1) 


{1,1} 


0.4687500 






0.4687500[2] 




6 


7(2) 


{2,0} 


0.0468750 


0.1482318 












{0,3} 


0.1482318 


0.4687500 




0.0000000[1] 


0.5156250[4] 


7 


8(1) 


{1,2} 


0.4921875 






0.4921875[4] 




8 


10(2) 


{2,1} 


0.0937500 


0.1960922 












{0,4} 


0.1960922 


0.4101562 




0.0000000 


0.5039062 [6] 


12 


19(3) 


{4,0} 


7.3242187xl0" 4 


1.0863572xl0~ 2 


1.5742811xl0 -2 










{2,3} 


1.0863572xl0~ 2 


0.1611328 


0.2335036 










{0,6} 


1.5742811xl0" 2 


0.2335036 


0.3383789 


0.0000000 


0.0000000 0.5002441[13] 



civ 



exact 

[i] 



$* BM . Since <f>f BM cx {z x - z c )(z 2 - z c )(z 3 - z c ) 
[see Eq. ||3J)], this result agrees again with the find- 
ings of Refs. [13, S|- For L = 5, the single non- 



spurious state is an excited one, $j^ act cx v 3 1^2- 

For L — 6, the ground-state is found to be <]>™ act oc 

-160$R BM /9 + Qi/4= (z 1 -^ 2 ) 2 (z 1 -z 3 ) 2 (z 2 -z 3 ) 2 ,i.e., 
the bosonic Laughlin function for v = 1/2 is equivalent to 
an RBM state that incorporates vibrational correlations. 
For L > N(N-1) (i.e., v < 1/2), the EXD yrast energies 
equal zero, and with increasing L the degeneracy of the 
zero-energy states for a given L increases. It is important 
that this nontrivial behavior [35| is reproduced faithfully 
by the present method (see Table Jj). 

Three electrons — Although unrecognized, the solution 
of the problem of three spin-polarized electrons in the 
LLL using molecular trial functions has been presented 



RBM/ 



Ref. 



33j. Indeed, the wave functions in Jacobi coor- 
dinates in Eq. (18) of Ref. [33| are precisely of the form 
QsPj as can be checked after transforming back to 
cartesian coordinates [36| . It is noteworthy that Laugh- 
lin did not present molecular trial functions for electrons 
with N > 3, or for bosons for any N. This is done in 
the present paper. We further note that the well-known 
Jastrow-Laughlin ansatz Y\i< i< j<N(zi — Zj) 2p+1 intro- 
duced in Ref. 0] is described as a quantum- fluid state Q • 
The quantum-liquid-picture interpretation of this ansatz 
is inconsistent with the RVM functions of Eq. (Q]) . 

Four electrons — For N — 4 spin-polarized electrons, 
one needs to consider rovibrational states [see Eq. (JTJ] for 
two distinct molecular configurations, i.e., $^ E |f (0, 4)Q™ 
and <Pf^{l,3)Ql l . Vibrational excitations with A > 2 
must also be considered. In this case the molecular basis 



states are not always orthogonal, and the Gram-Schmidt 
orthogonalization is implemented. Table |TT] summarizes 
the molecular description in the start of the LLL spec- 
trum (6 < L < 15 and L = 18). 

We note that in several cases the nonspurious states 
are given by a single trial state as defined in Eq. (JlJ. 
Indeed for L — 9 the yrast state is a pure REM state, 
i.e., <3?^ EM (1, 3). For L = 11 the single nonspurious state 
is the first excited state in the full spectrum, coinciding 
with the molecular vibration <I>f' EM (l, 3)Q 2 - 

Of particular interest is the L = 18 case; it corresponds 
to the celebrated v = 1/3 fractional filling, which is con- 
sidered Q as the prototype of quantum liquid states. 
However, in this case we found (see Table [TTJ) that the 
exact nonspurious solutions are linear superpositions of 
seven molecular states involving dipole (A = 2) and 
octupolc (A = 3) vibrations over both the (0,4) and 
(1,3) configurations. Focusing on the yrast state with 
L = 18, we found that its largest component is the 
pure $f 8 EM (0,4) REM state with a 0.9294 overlap with 
the EXD solution; the contibutions of the remaining six 
states are much smaller, but they bring the overlap to 
precisely unity. Unlike the v — 1/2 case of bosons, we 
stress that the fermionic Jastrow-Laughlin functions at 
all i/'s exhibit less-than-unity overlaps [3. \vi\ . 

Of great interest also is the L = 30 [v = 1/5) 
case, which in the composite-fermion picture was found 
to be susceptible to a competition 16[ between crys- 
talline and liquid orders. However, we found that the 
exact nonspurious states for L = 30 are actually lin- 
ear superpositions of the following 19 [= D Tl (L 
30)] RVM functions: 



*§S2(o,4)Q 



l 2 ~ 2k , with k 
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TABLE II: Spectra of four spin-polarized electrons. Second column: Dimensions of the EXD and the nonspurious TI (in 
parenthesis) spaces. Last three columns: Total energy eigenvalues [in units of e 2 /(«;Zs)] from the diagonalization of the 
Coulomb interaction e 2 /(Knj) in the TI subspace spanned by the trial functions "1 , ^™(0,4)Q™ and 3)Q™. Third 

to sixth columns: the molecular configurations (711,712) and the quantum numbers k, A and m are indicated within brackets. 
There is no nonspurious state with L = 7. An integer in square brackets indicates the energy ordering of each nonspurious 
state in the full spectrum, obtained by considering the spurious center-of-mass excitations (see also caption of Table [TJ> . Eight 
decimal digits are displayed, but the total energies agree with the corresponding EXD results within machine precision. 



L 


D EXD (£> TI ) 


[(ni,n 2 ){fc, A 


m}] 


Total energy eigenvalues (TI) 




6 


1(1) 


[(0,4){0,A,0}] 




2.22725097(1] 






8 


2(1) 


[(0,4){0,2,1}] 




2.09240211(1] 






9 


3(1) 


[(1,3){1,A,0}] 




1.93480798(1] 






10 


5(2) 


[(0,4){1,A,0}] 


I(0,4){0,2,2}] 


1.78508849(1] 


1.97809256(3] 




11 


6(1) 


[(1,3){1,2,1}] 




1.86157215(2] 






12 


9(3) 


[(0,4){1,2,1}] 


[(0,4) {0,2,3}] [(1,3){2,A,0}] 


1.68518201(1] 


1.76757420(2] 


1.88068652(5] 


13 


11(2) 


[(1,3){1,2,2}] 


t(0,4){l,3,l}] 


1.64156849(1] 


1.79962234(5] 




14 


15(4) 


[(0,4){2,A,0}] 


[(0,4){1,2,2}] [(0,4){0,2,4}] 


1.50065835(1] 


1.63572496(2] 


1.72910626(5] 






[(1,3){2,2,1}] 




1.79894008(8] 






15 


18(3) 


[(1,3){3,A,0}] 


[(1,3){2,3,1}] [(1,3){1,3,2}] 


1.52704695(2] 


1.62342533(3] 


1.74810279(8] 


18 


34(7) 


[(0,4){3,A,0}] 


[(0,4){2,2,2}] [(0,4){1,2,4}] 


1.30572905(1] 


1.41507954(2] 


1.43427543(4] 






[(0,4){0,2,6}] 


[(1,3){4,A,0}] Kl,3){2,2,3}] 


1.50366728(8] 


1.56527615(11] 


1.63564655(15] 






[(1,3){3,3,1}] 




1.68994048(20] 







.12-3^/2 



with k 



2,4,6; 

-k 



0,1,2,3,4,5,6; $£*$(1,3)Q* 

^ k (0A)Ql- 4k/ \ with k = 0,3; and *£^(1,3)Q§- , 
with k = 2,3,4,5,6,7,8. Diagonalization of the Coulomb 
interaction in the above TI subspace yielded an en- 
ergy 0.25084902 e 2 /(^s) per electron for the yrast 
state; this value agrees again, within machine precision, 
with the EXD result. The most sophisticated variants 
of the composite-fermion theory (including composite- 
fermion diagonalization (CFD), composite-fermion crys- 
tal (CFC), and mixed liquid-CFC states [II], M, ^ ) fall 
short in this respect. Indeed the following higher energies 
were found 0: 0.250863(6) (CFD), 0.25094(4) (mixed), 
0.25101(4) (CFC). The CFD basis is not translationally 
invariant. Consequently, to achieve machine-precision ac- 
curacy, the CFD will have to be performed in the larger 
space of dimension D EXD (L = 30) = 169. 

Conclusions. — The many-body Hilbert space corre- 
sponding to the translationally invariant part of the LLL 
spectra of small systems (whether fermions or bosons, 
and for both low and high angular momenta) is spanned 
by the rovibrational molecular trial functions introduced 
in Eq. fl}. The yrast and excited states for both short- 
and long-range interactions can always be expressed as 
linear superpositions of these trial functions. Thus the 
nature of strong correlations in the lowest Landau level 
reflects the spontaneous emergence of intrinsic point- 
group symmetries associated with rotations and vibra- 
tions of molecules of localized particles arranged in con- 
centric polygonal-ring configurations. We stress, particu- 
larly, the validity and numerical superiority of the present 
molecular theory for low angular momenta, where fun- 
damentally different " quantum- liquid" physical pictures 
(based on Laughlin, composite-fermion, and Pfaffian 



functions) have been assumed Q, [ll], [26[ to apply. 



[1 

[2 

[3. 
[4 
[5 

ir 

[8 

[9 

[10 

[11 

[12 

[13 

[14 

[15; 
[16 

[17 

[18 

[19 
[20 

[21 

[22 



D.C. Tsui, H.L. Stormer, and A.C. Gossard, Phys. Rev. 
Lett. 48, 1559 (1982). 

R.B. Laughlin, Phys. Rev. Lett. 50, 1395 (1983); Rev. 
Mod. Phys. 71, 863 (1999), and references therein. 

F. D.M. Haldane, Phys. Rev. Lett. 51, 605 (1983). 

B. I. Halperin, Phys. Rev. Lett. 52, 1583 (1984). 

S.A. Trugman and S. Kivelson, Phys. Rev. B 31, 5280 
(1985). 

J.K. Jain, Phys. Rev. Lett. 63, 199 (1989). 

G. Moore and N. Read, Nucl. Phys. B 360, 362 (1991). 
M. Stone et al, Phys. Rev. B 45, 14 156 (1992). 

J.H. Oaknin et al., Phys. Rev. Lett. 74, 5120 (1995). 
J.J. Palacios and A.H. MacDonald, Phys. Rev. Lett. 76, 
118 (1996). 

J.K. Jain, Composite Fermions (Cambridge University 
Press, 2007), and references therein. 

C. Yannouleas and U. Landman, Phys. Rev. B 66, 115315 

(2002) . 

C. Yannouleas and U. Landman, Phys. Rev. B 68, 035326 

(2003) ; ibid. 70, 235319 (2004). 

C. Yannouleas and U. Landman, Rep. Prog. Phys. 70, 
2067 (2007). 

C.C. Chang et al, Phys. Rev. Lett. 94, 016809 (2005). 
C.C. Chang et al, Phys. Rev. B 73, 155323 (2006). 
B.R. Mottelson, Phys. Rev. Lett. 83, 2695 (1999). 
K. Wilkin and J.M.F. Gunn, Phys. Rev. Lett. 84, 6 

(2000) . 

N.R. Cooper et al, Phys. Rev. Lett. 87, 120405 (2001). 
Th. Papenbrock and G.F. Bertsch, Phys. Rev. A 63, 
023616 (2001). 

M. Ueda and T. Nakajima, Phys. Rev. A 64, 063609 

(2001) . 

M. Popp et al, Phys. Rev. A 70, 053612 (2004). 



5 



[23] C.C. Chang et ai, Phys. Rev. A 72, 013611 (2005). 
[24] N. Barberan et ai, Phys. Rev. A 73, 063623 (2006). 
[25] L.O. Baksma ty et al, Phys. R ev. A 75, 023620 (2007). 
[26] N.R. Cooper. larXiv:0810.4398V l. 

[27] N.D. Gemelke, Ph.D. Dissertation (Stanford University, 
2007); Steven Chu, Book of Abstracts, ICAP 2008, p. 5. 

[28] An yrast state has the lowest enegy at a given angular 
momentum. 

[29] G.S. Jeon et ai, Eur. Phys. J. B 55, 271 (2007). 
[30] C. Yannouleas (unpublished) . 

[31] For bosons, Ref. [201 ] studied some aspects of the Q\ ex- 



citations in the range < L < N. 

[32] S. Viefers, J. Phys.: Condens. Matter 20, 123202 (2008). 

[33] R.B. Laughlin, Phys. Rev. B 27, 3383 (1983). 

[34] R.A. Smith and N.K. Wilkin, Phys. Rev. A 62, 
061602(R) (2000). 

[35] X.-G. Wen, Quantum Field Theory of Many-Body Sys- 
tems (Oxford University Press, 2004). 

[36] Use of an algebraic language like MATHEMATICA is 
necessary. 



